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Abstract 

We investigate analytically and numerically the spatial structure of the non-equilibrium 
stationary states (NESS) of a point particle moving in a two dimensional periodic Lorentz 
gas (Sinai Billiard). The particle is subject to a constant external electric field E as well 
as a Gaussian thermostat which keeps the speed |v| constant. We show that despite the 
singular nature of the SRB measure its projections on the space coordinates are absolutely 
continuous. We further show that these projections satisfy linear response laws for small E. 
Some of them are computed numerically. We compare these results with those obtained from 
simple models in which the collisions with the obstacles are replaced by random collisions. 
Similarities and differences are noted. 



1 Introduction 

In this paper we continue our study of nonequilibrium stationary states (NESS) maintained 
by a Gaussian thermostat [U HI EJ [9] . Theoretical analysis and computer simulations show 
that the NESS obtained from such artificial model dynamics can give useful information on 
real systems maintained in NESS by coupling with heat baths [8]. 

Here we focus on the Moran-Hoover (MH) model of a single particle in a periodic billiard 
moving under the influence of an electric field E and a Gaussian thermostat that keeps the 
kinetic energy constant [12]. The equations of motion are: 



x = p 

p = E - a(p)p + F obs (x) (i) 

n ( n \ - (P-E) 



where x is the position, p the momentum of the particle with unit mass, and F D b s (x) repre- 
sents elastic collision with the obstacles. It is clear from eq.([T|) and the fact that collisions 
with the obstacles do not change |p| that 4(p ■ p) = 0. We shall therefore set |p| = 1 from 
now on. 
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The particle moves on a 2-dimensional torus whose side can be chosen to be 1. An 
arrangement of the obstacles used for all the simulations presented in this paper is shown in 
Figure [TJ The two obstacles have radii r\ = 0.2 and r2 = 0.4. This is also the arrangement 
used in our previous works [TJ [2] and is chosen to have a finite horizon, i.e. there is an upper 
bound for the time between successive collisions of the particle with the obstacles. Moreover 
we take E to be along the horizontal x axis, i.e. E = (E, 0). The analytical results apply to 
general geometries with finite horizons. 



E 




Figure 1: Typical obstacles placement. 



The set of states where the particle collides with a given obstacle can be parametrized by 
two angles: i? £ [0, 2ir] the angle on the obstacle between the collision point and the positive 
x direction, and ip the angle between the particle velocity and the outgoing normal to the 
obstacle at the collision point. To obtain a complete coordinate system for the collision states 
we define the coordinate 9 = f? for obstacle 1 (see Figure [1]) and 9 = i? + 2tt for obstacle 2 
so that 9 € [0,47r]. In these coordinates, the elastic collision is simply represented by the 
map C(9,tp) = (9,ir — ip), where ip € [7r/2,37r/2] before collision and ip € [— 7r/2,7r/2] after 
collision. We will call M = [0, 4tt] x [— 7r/2,7r/2] the set of possible pairs (9,ip) representing 
the position of the particle just after a collision. Ai corresponds to a Poincare section of the 
flow. See Figure [2] for a depiction of 9 and ip. 

Since |p| is constant, the trajectory of the particle can be represented by its position x(i) 
and the angle of its momentum with the horizontal axis <p(t). The motion of the particle be- 
tween two collisions can be exactly integrated. Moreover one can construct the map Se(9, ip) 
mapping the position and momentum of the particle just after a collision to its position and 
momentum just before the next collision, which may be with the same or a different obstacle. 

In this way we can represent the dynamics in discrete time as the iteration of the map 
Te '■ M — > A4 between successive collisions given by Te = C o Se- Observe that this map is 
not continuous (grazing collisions) and that, for E small, the dynamics is a perturbation of 
the free billiard dynamics [4]. 

In our previous works we were primarily concerned with the SRB distribution associated 
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Figure 2: Elastic collision. 



with Te- Let //o be the measure on Ai given by 

im = caa(il>)x(9)d6dil>/Z (2) 

where Z = 47r(ri +r2) is a normalization constant and x($) = r i f° r G [0, 27r] and x(#) = ?~2 
otherwise. Observe that /j,q is invariant under To. The SRB distribution /iE(dO, dip) is defined 
as the weak limit of /j,q under the dynamics Te, i.e. 

i x E = w-lim Tg/u . (3) 

The measure he, when it exists and is unique, represents the natural non equilibrium steady 
state (NESS) for the system [13] . Clearly /j,q is the SRB measure of Tq. 

From the SRB measure ^E(dO, dip) on A4 for the collision map one can build the SRB mea- 
sure mE(dx, dip) on M = Q x [0, 2ir] for the flow generated by eq.([T]). Here Q is T\ obstacles. 
This can be represented as: 

m E (A) = ^- / / I A (X?(6,iP),$F(6,iP))dtv E (d9,diP) (4) 

te Jm Jo 

where Ia is the indicator function of the set A C M, (Xf,$>f) is the flow generated by 
eq.(pQ) and te(0,iP) is the time till the next collision when starting at (0,ip) € M with 
te = j M te(0, ip) ^E{d0, dip) denoting the mean free time. 

In [4j it was proved that, for small fields E, \E\ < Eq, the above model has a unique NESS 
described by an SRB measure \ie which is singular with respect to the Liouville measure with 
Hausdorff dimension given by the standard Kaplan- Yorke formula [llj . 

The current j(E) in this NESS is given by 

j(E) = m E (v) 

where v = (cos (</>), sin (0)) is the velocity of the particle. This current was shown in [1] to 
be given by the Kawasaki formula cf. [6]. In the limit E — > the Kawasaki formula reduces 
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to the Green-Kubo formula for the conductivity n which satisfies the Einstein relation [HE]. 
An investigation of the current as a function of the field was carried out in [1] . It was argued 
there that the current is not a C 1 function of the field E even close to E = 0. The results of 
[3] were generalized in [H [10] to systems where the collision rule or the free flow dynamics is 
perturbed. 

In none of the above works was the spatial dependence of the singular (with respect to 
Lebesgue) measure m E (dx,d<j)) studied. This is what we do in this note. We will describe 
analytical results and numerical studies of the spacial and angular dependence of the NESS 
m E (dx,d(f>) when projected on <f> G [0, 27r] or on x G Q and related quantities like the local 
average velocity. 

The outline of the rest of the paper is as follows. In section [2] we introduce the local 
density, local average velocity and angular distribution derived from m E . We find their 
dependence on position and field strength. We also show there computer generated pictures 
of the flow and compare them with the predictions of Green-Kubo formulas at small E. In 
section [3] we introduce and analyze two simple models in which the deterministic collisions 
with fixed obstacles are replaced by random collisions whose times form a Poisson process and 
compare their properties with those of the deterministic model. The appendixes are devoted 
to analytical justification of the claims in section [2j A paper describing results for the case 
where the system consists not just of one but of a large number of particles is in preparation. 



2 Local Structure of the SRB measure 

We define and study the several projections of the SRB measure rriE introduced in the previous 
section. For clarity of exposition we delay derivations and justifications to the Appendices. 



2.1 Local Density and Average Velocity 

Two interesting quantities to study are the local density and local average velocity. More 
precisely, we define the projected measures on the position x as: 



5 E (A) = / m E (dx,d<f>) (5) 

JAx[0,2n) 

for any A C Q. This clearly defines a probability measure 5 E (dx.) on Q. Using eq.Q, and 
defining 

J${0,tI>) = —1 I A {X?(B,il>)) dt (6) 

T E Jo 

where I a is the indicator function of the set A, we can represent 5 E (A) as the integral of a 
piecewise smooth function with respect to u. E (dO,d(f)): 

5 E (A)= [ j%{e^)n E (de,dip) (7) 

J M 

Observe that J^(B, is the relative amount of time the trajectory starting from (6, i/j) spends 
in the set A, that is the amount of time divided by the mean free time f E . 
We also define the vector measure for the local average velocity 



v E {A)= I (cos(0),sin(^))m E (dx,d0) (8) 

JAx[0,2n) 



1 



Also this measure can be written as the integral of a piecewise smooth function with respect 
to fiE(d9,dcj)), see eq. ([52|) below for details. 

In Appendix IA1 we show that, for \E\ < Eq, the integrals in eq.© and eq.(|8]) define 
absolutely continuous measures. That is <5e(dx) = n^(x)dx and z/£:(dx) = 7?,£:(x)ve(x) dx. 
We call ng(x) the local density and ve(x) the local average velocity at x. We show that both 
are continuous functions of both x and E with n-o(x) =const= (Area(Q)) 1 and vo(x) = 0. 

To visualize the above numerically, we divided the torus of Figure Q] in a grid of 50 x 50 
cells and computed the time average of the velocity of the particle when it crosses a cell. The 
results are shown in Figure El We also computed the local density on the same grid. 




Figure 3: Average local velocity v#(x) for E = 0.1. 

We now show that the local density ue(x) and the local average velocity vg(x) are linear 
in the field E when E — > 0, that is 

n E (x) = n (x) + d(x)E + o(E) (9) 
Vij(x) = k(x)E + o(E) (10) 

where d(x) and k(x) can be computed via Green-Kubo-type formulas as follows. 

Consider the family of all velocity vectors originating at the point x (at which we are 
computing the density or average velocity) ; they make a one parameter family of phase states 
W x = {(x, (f>) | < (j) < 2ir}. Let p x be the probability measure on Wx that has a uniform 
distribution over </> G [0, 2ir]. We can map VF X to the collision space A4 by taking every point 
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(x, (j)) G W x to its first collision with dQ in the past, under the field-free dynamics. The image 
of W x will then be a collection Wo of curves in Jvi on which we get an induced probability 
measure pq. Pulling this measure further back (into the past) we get a sequence of probability 
measures p n = T -n (/5o), each sitting on a collection W n = T -n (Wo) of curves in M. With 
this definition we get that 



d(x) 



Po(A ,x) + ^Pn(A 



n=l 



fill 



where c = Area(Q) -1 = no(x), Ao = ro(#, VO cos (^ + iji) is the x-distance form the collision 
point {0, tp) G A4 to the next collision point and Ao jX is the x-distance from the collision 
point (6, ip) G .M to the point x. Observe that po is supported on points whose trajectory 
passes through x before colliding again. The above series converges exponentially, because 
the measures p n converge exponentially fast to the measure po on A4 (see Theorem 7.31 in 
0) and ^o(Ao) = 0. 

Consider now the two signed measures p c ' x and /3 S,X on W x that have densities cos 4> an d 
sin<^>, respectively, with respect to p x . As before, we can map p c ' x and p s ' x on the collision 
space M. and obtain signed measures Pq and p s Q on Wo, respectively. We also denote their 
images by p% s = T ~ n (/5Q S ) on W n for nCZ, 

oo 

k(x) = i (Pn(Ao),^(A )). (12) 



The terms in the series in eq. (|12j) converge to zero as n — > ±oo exponentially fast, because 
the measures p% s converge to the zero measure; this again follows from Theorem 7.31 in [5]. 

We note that the perturbation of the density n^(x) and of the local average velocity 
ve(x) are linear in E, to the leading order, and the factor of E is given, in both cases, by an 
infinite sum of correlations, i.e. the right hand sides of eq. (|1 lj) and eq. (|12p . 

We computed numerically the coefficients d and k in eq.([9|) and eq. tllOp to compare their 
predictions with the simulation results shown in Figure [3j We truncated the infinite sums 
in eqs. (|llll2p to \n\ < 15 since we saw no visible difference arise from taking more terms 
into consideration. Let l£ = {(x,y) G Q} be the vertical cross section placed at horizon- 
tal coordinate x and l~ = {(x,y) G Q} be the horizontal cross section placed at vertical 
coordinate y. Finally let e x = (1,0) and e y = (0,1) be the unit vectors in the horizontal 
and vertical direction respectively. Figure H] shows a comparison of the horizontal component 
(ve(x) • e x ) of Vg(x) along Zq" 41 with the prediction of eq. (fl0|) . In the same way, Figure[5] 
shows a comparison of the vertical component (vg(x) • e y ) of vg(x) along /q~ 41 again with 
the prediction of eq. (110|) . In both figures the pluses represent the results of direct simulation 
while the crosses are obtained using the Green-Kubo formula eq . (|12|) . 

The comparison of ng(x) with eq.([9]) is more difficult. Calling n^(x) = (n#(x) — 
71_e(x))/2 and n e E (x) = (n^(x) + n_s(x))/2 — no(x), we have that n° E {x) satisfies the 
same linear response formula eq.© of n^(x) with the same coefficient d(x) but we expect 
the remainder to be smaller. This is relevant in the present case since n E (x) and n E (x) ap- 
pear to be of comparable magnitude. We observe that, due to the symmetry of the problem, 
ng(l — x, y) = u-e(x, y) so that n E (x, y) = (ue{x, y) — n^(l — x, y))/2. Figure [6] compares 
n E (x) along with eq.([9|). Again the pluses represents direct simulation while the crosses 
are obtained using the Green-Kubo formula eq. (fTTI) . 

More generally, given a probability measure p(dx,d(p) = l(x,cf)) dxdcf) absolutely continu- 
ous with respect to the Lebesgue measure on Q let pf{dx.,d4>) = lf{x.,4>)dx.d4> be its time 
evolution with respect to the dynamics generated by eq.([T]). In a similar way as above, we 
can then define: 
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Figure 4: The x component of the average local velocity vg(x) for E = 0.1 and x = 0.41. 



«f(x) 
«f(x)vf(x) 



The density nf (x) clearly satisfies a conservation law: 



(cos(</>),sin(</>))/f (x, 



| / nf (a) dx = - [ nf (x)(vf (x) ■ n(x)) dtr(x) 



(13) 



(14) 



where ^4 is a subset of Q with smooth enough boundary, n(x) is the unit outward normal to 
dA at x and <r(x) is the length element on dA. 

Taking the limit t — > oo and assuming that lim^oo nf(x) = n^(x) and lim^oo vf (x) = 
Ve(x) we obtain 



/ n E (x)(v s (x) • n(x))dcr(x) = 
JdA 



(15) 



The above assumption is not trivial. It is easy to show that, if lim^oo nf{x) exists, it 
has to equal ng(x). On the other hand, we do not have a proof for the existence of such 
a limit. A similar argument holds for wf (x). A complete justification of eq. (|15j) will thus 
require further work but we certainly expect it to be true. 

Nonetheless we can test the validity of eq. (|15j) numerically. Due to the symmetry of Q 
we have that the average current j(E) = (j(E),0). Moreover, since the collision are elastic, 
v#(x) is tangent to dQ for x G dQ. It follows from this that 



/ n E (x)(v B (x) • e x ) dy = j(E) 

J l x 
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Figure 5: The y component of the average local velocity vg(x) for E = 0.1 and y = 0.41. 

/ n s (x)(v s (x) • e y ) dx = 

Jly 

independently on the value of x or y. Both these equations are very well verified. 

2.2 Angular Distribution 

We now look at the projection of mg on the angle eft. We can define the projected measure 
rj(d(j)) by setting, for a measurable set A C [0, 2ir], 

rj E (A) = / I A (x, 4>)m E {dx., #) 
Jm 

where I a is the indicator function of the set A. Again we can write i]e(A) as an integral on 
the SRB measure ^E(d6,dip) as follows. Define the function: 

Ja(W) = — I A (*f(9,il>)) dt (16) 

T E Jo 

Then we have that 

"He (A) = [ J%(9, ^) v E (dd, d<f>) = n E (jf) 
Jm 

Using the argument in Appendix |A] we can show that, for \E\ < Eq, rj E is absolutely contin- 
uous with respect to d(f>, i.e. that rjE(d4>) = hE(<j>)d<f) where hE{4>) is a continuous function 
of both <j) and E with ho((j)) =const= l/2n, since the invariant measure uiq is uniform on 
Q x [0,2vr]. 

We computed h,E{<p) numerically for E = 0.1. The result is shown in Figure [3 A striking 
characteristic of this distribution is the dip around = 0. This is somewhat unexpected since 
the effect of the field E is to push the velocity of the particle to align with the positive x 
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Figure 6: The symmetrized local density n° E (x.) for E = 0.1 and x = 0.41. 

direction so that one would expect a maximum at (ft = rather than a local minimum (see 
also section [3] for a comparison with the stochastic models). 

To understand this better we consider, for a given (ft, all points (9, ift) S M. that produce 
the outgoing velocity vector (cos (ft, sin (ft), i.e., we consider 



VA = {(p,iJ,)eM:il> + 



(mod 2tt)}. 



Now A4 is foliated by the lines {V^,}, < (ft < 2ir. Let /Xq denote the conditional measure 
induced by hq on the line V^. If we use 9 as the (only) coordinate on Vd>, then 

df4 = Z- 1 cos(<A - 6)x(0) d6 
where Zm is the normalizing factor 



Z^= / cos(0-0)x(0)d0 = 2(ri+r 2 ). 

/cos(</>-6»)>0 



(17) 



We remind the reader that x(^) = ri and < 9 < 2tt on the first obstacle and x(^) = r 2 an d 
2-7T < < 4-7T on the second. 

Now we consider the conditional distribution of the free flight time function To on each 
line Va). It turns out that its first moment is constant, i.e., /^o( r o) = ^0 f° r a H <^'s, where 

— / / o( r ) is the total (unconditional) mean free time. In other words, the deterministic 
collision process is isotropic, on average. This seems to be a novel result in the studies of 
billiards and we prove it in Appendix We now argue that the observed dip near (ft = 0, 
for small E can be traced to the second moment, ^q(tq), which is not constant and which 
for our obstacles indeed has a local minimum at (ft = 0. 

We will show that the density 1ie{4>) satisfies 



h E {(ft) 



1 

2^ 



+ a{(ft)E + o{E) 



(18) 
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Figure 7: Angular distribution for E = 0.1. 



where a(4>) is given by a Green-Kubo formula 



a{4>) 



7 1 00 

YW Q ^ 4(ro-(A oT ")). 



(19) 



Recalling that Z = 4ir(ri + r2) is the normalization of fio, see text after eq.©, and Z^ = 
2(n + r2) is independent of cj>, see eq. (fT71) . we have that Z^/Z = l/2n. Again we see that 
the fluctuations of the density He are linear in E, to the leading order, and the factor of E 
is given by an infinite sum of correlations. The latter converges exponentially fast according 
to general results (Theorem 7.31 in [5]). 

Usually its central term (n = 0) is the most significant, and it is given by 



cos q 



(20) 



The central term explicitly involves the second moment of to restricted to V^. Even 
though cose/) has a maximum at eft = 0, it may be more than counterbalanced by a dip that 
the second moment A*q ( r o ) ^ as nea r = 0. This is exactly what happens in our model shown 
in Figure [TJ 

To check numerically the above results we proceed like in the case of n^(x) in Figure 
[3 We introduce the odd part of the angular distribution h E (cp) = (hE(4>) — /i-b(</>))/2 and 
observe that it satisfies the linear response equation 



h E {ck) = a{ck)+o{E) 



{21) 



with a(4>) still given by eq. (|19p . Again we expect the reminder to be smaller. Finally, due to 
the symmetry of our system, we have that h° E ((f)) = (/i£;((/>) — /i^(0 + tt))/2. 

Figure [8] presents the plot of eq. (f2T|) and the numerically computed plot of h E (4>) for 
E = 0.1. The crosses represent the numerically computed value of h E (4>). The pluses come 
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from the central term of eq. (|19p . Already at this level the dip is clearly visible and the 
agreement is pretty good. Finally the boxes represent eq . (|19|) truncated at |n| = 20. We 
have computed eq. (|19j) truncating the sum up to \n\ = 100 but no significant difference from 
|n| < 20 can be observed. This is clearly consistent with a fast convergence in the sum in 



h° 




Figure 8: Comparison between the angular distribution for E — 0.1 and eq.f l59p ; see details after 
eq.flSB. 



Our analysis indicates that the dip at <j) = appears to be an artifact of the geometry of 
the scatterers chosen for our deterministic model. 



3 Random Collision Models 

In [2], we introduced a simplified version of the MH model by replacing the collisions with 
the fixed obstacles with a Poisson random collision process. More precisely we assume that, 
in every time interval dt the particle has a probability equal to A|p|di (with |p| in this case 
fixed to be 1) to undergo a collision. Between collisions the particle moves according to eq.([T|) 
without -Fobs- When a collision happens we consider two collision rules: 

I the velocity of the particle after the collision is in direction </> € [0, 2ir] with probability 
density d<f>/2'K\ or 

II an angle 7] £ [— tt/2, it/2] is chosen at random with probability proportional to cos(r])dr] 
and the direction of the velocity is changed according to an elastic collision rules for a 
particle colliding with an obstacle with outgoing velocity forming an angle r\ with the 
normal to the obstacle. 

We call the models with the above collision rules Model I and Model II. 

We can think of Model II as representing a situation in which we have N scatterers 
with diameter e randomly placed in T and we consider the (Boltzmann-Grad) limit in which 
N — > oo, e — > 0, such that iVe 2 — > while iVe — > A" 1 , the mean free path [7]. 
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Let f a (E,x,(j),t) be the probability density at time t of finding the particle at x with 
momentum p = (cos (ft, sin (ft). Here a =1,11 indicates Model I or Model II respectively. This 
density satisfies the equation: 

d t f a (E,x,(ft,t) - pd^f a (E,x,(ft,t) - Ed e (sin 9 f a (E, x, (ft, t)) = 

A ^ p a (v)fa(E, X, (/) + 2r/ + 7T, t) dl] — f a (E, X, (ft, f)^ (22) 

where, i G M+, xET, the unit torus, E = (E,0) is in the horizontal direction and A is the 
collision rate. Moreover we have pi{rj) = it" 1 for Model I and pii(rj) = cos 77/2 for Model II. 

It follows from eq. ([22|) that, when the distribution at time 0, f a (E, x, (ft, 0) does not depend 
on x, the density f a (E,x,(ft,t) will also not depend on x for every t > 0. Even if the initial 
state does depend on x, it is easy to show [3] that as t — )• 00 the system will approach a 
stationary density f a (E,(ft) which will satisfy the equation: 

7T 

- jd^ (sin 4>f a (E, (ft)) = J^p a ( V )f a (E,4> + 2ri + iT)dri-f a (E,(ft). (23) 

From now on we will set A = 1 since the stationary f a depends only on E/X. 

We can try to solve this equation as a power series in E. Since E is a singular perturbation 
the series will not be convergent for any non zero value of E. However, we expect that it will 
be an asymptotic series and accurate for small \E\. Writing 

00 

UE,(ft) = Y. El ti ] ^) (24) 
yields a hierarchy of equations for i = 0, 1, 2 . . . : 



- (sin^- 1 )^)) = / \ Pair)) (eft + 2 V + n)d V - (25) 

with fa ^ = 0. The equation for i = is easily solved and gives, as the unique solution, 
fa^ = 1, since we require f f a (E, (ft) d(ft = 2-k. To solve the higher order equations we write 



fi i] {4>)= £ /«(n)cos(n^) 

n=— 00 

where we used the symmetry with respect to the direction orthogonal to the field to eliminate 
the terms in sin(n6>) and clearly fa\n) = fa\—n). In this way, for n ^ 0, eq. ([25j) becomes 

= 2 (l-p a (n)) (ft 1] (n - 1) - /(«-D(n + 1)) (26) 

with pi(n) = for Model I and pn(n) = l/4n 2 for Model II. Finally /« (n) = (5 n> o, again due 
to the normalization condition. This yields 



0) = 1 + E cos{(ft) + cos(2<£) + . . . 

3 45 

/„(£;,</)) = l + - J Bcos(^) + — J E 2 cos(2«/)) + ... (27) 

4 64 

We can compare the above results with numerical simulation of the stochastic processes 
generating eq. (j23|) . We set E = 0.2 and run both processes for 10 s collisions. The results 
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Figure 9: Comparison between numerical simulations of the stochastic process and the power series 
expansion eq. (j27J) for E = 0.2. See explanation after eq. (j2"Tj) . 



are plotted in Figure [9j The crosses refer to Model I while the pluses refer to Model II. 
Superimposed are the graph obtained from eq. (|27|) . As one can see, the fit is very good. This 
is in agreement with our expectation that the series in E is an asymptotic one. In the case 
of model I this can be rigorously justified, see eq. fpUj) below. 

We note that, for both Models, the power series for / a (0.2,0) has a global maximum for 
= 0. Since this is not true for the angular distribution of the deterministic MH model, see 
section [2T2l we will investigate the behavior of f(E,ip) near = more closely. 

3.1 Model I 

Eq. (p3|) can be written as: 



-Ed 4> {sin<j>f l {E,<t>)) = l-fi{EA) (28) 

where we normalize f\ as j fi(E, 0) d<f> = 2ir. 

Eq. ([28|) can be solved by introducing the function 



h(E, 



1 - cos <A 2s 

1 + COS J 

h(E,<t 



which is a solution of the differential equation d^h(E, 0) = g , and defining 



sinffl 



Substituting in eq. ([28j) we obtain 
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d4>g{E,if>) = -h{E,4> + -n) 

Observe that h(E, 4>) has a non integrable singularity at <p = tt so that, for fi(E, <f>) to be 
integrable we need g(E, tt) = 0. We can thus represent the solution as: 

h(E, 0) = r h(E, rj + Ti)^ (29) 



£/sin(0) 

We list below some properties of fi(E,(f>) that will be useful in the following. We have two 
possible situations: 

E < 1 : In this case fi(E,(p) is continuous in (f> for every <j). Moreover it is easy to see that 
fi(E,(f)) is C°° for 4>^0,tt. For <f> = 0, tt, if E < l/n, fi(E,</>) is C n ~ x and d$fi(E,<f>) 
is Holder continuous of exponent a for < a < 1/E — n. 

E > 1 In this case fi{E,(j)) is still C°° everywhere but for cf) = 0,tt. At (j> = we have 
a singularity and fi(E,<f)) ~ 1 '- B_1 . More precisely the function fi(E,(f>) is 

Holder continuous of exponent < a < 1/E — 1/E' ', for every £" > i?. 



Starting form eq. (|29j) and integrating by part we obtain: 

/:(£?, 0) = 1-^^ ^ cos V h(E,r l + TT)d V = 
sm(<f>) 

= 1 + Ecos<j) + E k<yE ^) f cos2ri h(E,r> + tt) d-n = 

= l + £cos^ + g 2 cos(20) + E 2 kiyE ^} [ d (j ,[cos2r] sin r]]h(E,r] + tt) drj = 

sm(4>) 

N 

= Y; En fl l) ^)+ ENR N(E,<t>) (30) 
«=o 

The above expansion coincides with the one obtained in eqs. (|24H2"6]) . It is not difficult to see 
that \Rn(E, (f>)\ < KC N N\. Since it is clear from eq. (|29p that fi{E,<f>) is not analytic in E 
for small E, this inequality means that, as we discussed previously, the perturbative series 
for fi(E, (p) is at least asymptotic. Notwithstanding this, eq (l30p and the regularity properties 
of fi(E, 4>) tell us that, for E small, fi(E, 0) has a unique maximum at <f> = and a unique 
minimum at (j) = tt. 

3.2 Model II 

We can use the solution of Model I to get more analytical information on Model II. Proceeding 
as in eq. ([28j) we write the solution of eq. ([23j) as 



E sm <p 

and obtain the representation for g(E,<p): 

d^g(E, 0) = ~-^h(E, + tt)J cos 
from which, reasoning as in Model I, we get 



h(E, uj + tt) , „ . , . , 

\ ^ M E,CJ + 7T)d0J 31 

sm(cj + tt) 
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ME, 



1 



4E sin( 



-h{E, 



h(E,rj + 7r) 



cos 



00 — 7] 



f(E,uj + ir)dujdr] (32) 



for < (ft < tt. We can then set fn(E 
be written has as 



where 



00 



>) = fn(E,(j)). Observe that the above equation can 

fu(E, (ft) = r Q(<f>, oo)f n (E, w) duj (33) 

i r w 

h(E,(f>) / h(E, V + iT) 



for < (ft < tt and Q(- 
Moreover we have 



Q(4>,uo) d(ft 



AE sm(4>) 



,oo 



8111 



UJ — 7] 



drj 



oo). It is easy to see that Q(<p,uj) > for every 4>,uj. 



lim 

e->0 

1 



Q{(j),uj) d(f> = 
h(E,4>)h(E,(f) + tt) 



o 



sin 



co — (f) 



-Yin\h(E,e) / h(E,r] + 7r) 



sin 



LO — 7] 



dij 



+ 



— lim h(E, tt — e) 
4 e^O 7 



h(E,r} + ir) 



sm 



to — r] 



dr] 



1 

4.,, 



sm 



u) — 4> 



(34) 



where we have used that 4>) = h(E, cf> + tt) 1 and that, for e small, /i(-E, e) ~ e l / E while 
h(E, e + tt) ~ e" 1 /^. Proceeding in the same way for — 7r < < 0, we get 



Q{(f),u)) d(j) = - 



sm 



00 — (j) 

2~ 



1 



for every to. Finally, in the same way we got the regularity properties of fi(E,(p), we can 
see that, if E < 1 then Qfu is a Holder continuous function with Holder norm bounded by 
the L°° norm of f\\. This immediately implies, by the Ascoli-Arzela theorem, that Q is a 
compact linear operator on C . In this situation we can apply the Krein-Rutman theorem, see 
|15} I16j. and obtain that there is a unique function fu(E,(p) that satisfies eq, (|33p . Moreover 
fn(E,(f>) > for every <f). A similar argument tells us that, for E > 1, there is a unique 
solution of eq. (133p and it can be written as 

/„(£J,0) = |sin(0)| 1 -iz(E,^) 

with l{E,(f)) continuous in (ft and strictly positive. 

Observe that, for any integrable function f{(ft), we have 



cos 



00 — (ft 



f(uo + Tr)doo 



sgn(w — (ft) 



sm 



00 — (ft 



f(E, oo + Tr)doo 



2 



cos 



00 — (ft 



f(E,oo + Tr)doo 



COS 



00 — (ft 



f(E,U0 + TT)dG0 + f((ft) 



Thus the above integral is always at least C , while if f((ft) is C n it is C n+2 . This implies that 
fn(E,(ft) has the same regularity properties as a function of (ft as f\(E,(ft). In particular if 
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E < 1/3, fn(E, 4>) is C 2 and we can try to compute f{{{E, 0) explicitly. Observe that eq. (|23l) 
tells us that 



E sin 4>f{\(E, 4>) + E cos <j>f u (E, , 
Evaluating at <fi = we get 



cos 



oj — (j) 



fu(E,cu + 7T)(ko + fii(E,(j)) (35) 



fll(E, 0) = l _ — J* cos (|) f n (E, u + tt)^ 



(36) 



As expected, this equation loose meaning when E > 1 since fu(E,(fi) is no more continuous 
at = 0. We can now differentiate both side of eq.(|23p and obtain, after evaluating in <j) = 0, 

(1 - 2E)f{ 1 (E,4>) = ^J* sin (|) f n (E,uj + 7r)doj = 

for symmetry reasons. Again this equation make sense only if E < 1/2. Finally, differentiating 
once more, we get 



3Ef{[(E,0) - Ef u (E,0) = ^f cos (|) fn(E,u>)du - \fn(E,0) + f{[(E,0) 
Using eq. (|36|) we get 

fl[(E,0) = -^ T ^f u (E,0) (37) 

that is clearly negative for E < 1/3 so that we have that f(E,<fi) has a local maximum at 
<fi = 0. Observe that expanding this formula to third order in E we get a result in agreement 
with the expansion in eq. (!27p 



Acknowledgment 

The authors thank Eric Carlen for many insightful comments and discussions. The work 
of FB was supported in part by NSF grant 0604518. The work of NC was supported in 
part by NSF grant DMS-0969187. The work of JLL was supported in part by NSF grant 
DMR08021220 and by AFOSR grant AF-FA9550-07. The authors are also grateful to the 
Alabama supercomputer administration for computational resources. 



Appendices 

A Regularity of Projections of SRB measures 

SRB measures are characterized by absolutely continuous conditional distributions on unsta- 
ble manifolds, but generally they are singular. Singularity of a measure \x means that there 
is a subset Ai' C M. in the phase space M. such that fi(Ai') = 1 but the Lebesgue measure 
of M' is zero. In that case /i does not have a density on Ai. 

However in physics one rarely observes measures on the entire phase space; it is more 
common to observe distributions of some selected variables (e.g., positions or velocities of 
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selected particles). The distributions of those variables are obtained by projection of the 
relevant measure onto the corresponding variables. And the resulting distribution is often 
absolutely continuous, with a continuous density, despite the singularity of the measure in 
the whole of phase space. 

Similar smoothness results hold if, instead of projecting the measure onto certain vari- 
able^), we integrate some smooth functions with respect to all the other variables. For 
example, in eq.® we integrate cos</> and sin0 with respect to <f> and get an absolutely con- 
tinuous distribution with respect to x with a continuous density. 

We sketch a proof here that in our Moran-Hoover model the corresponding projections 
have continuous bounded densities. Let pe denote the SRB measure on the collision space 
with coordinates (9,ip) for a given value of E. Since we will consider only a given value of 
E, we will suppress the dependence on E in what follows. All estimates are uniform in E. 
Consider a projection of p onto a line transversal (not parallel or perpendicular) to stable 
and unstable manifolds and singularity manifolds. For simplicity, let p be projected onto the 
9 axis. 

The density of the projection can be computed as 

p{9) = lim <5-V(#0,5) 

<5->0 

where Rg t s = {{6' , tp) : < 0' < 9 + 5} is a rectangle in the collision space of size 5 in the 9 
direction. It is known that the SRB measure p satisfies 

p{r(6,ip) < 5) < C5 V5>0 (38) 

for some constant C > which is uniform for all small fields E, see [9]. Here r(9,ip) denotes 
the distance from (9,ip) € M to the nearer endpoint of the unstable manifold passing through 
(9, tp) (that manifold is a smooth curve which is divided by (9, ip) into two segments; so r(9, ip) 
denotes the length of the shorter one). 

Now for each unstable manifold W, the intersection W n Rg s is a segment of W that has 
length bounded above by C'6 for some constant C > 0. Hence its relative measure (within 
W) is of the same order as the measure of the segment of length 5 at an endpoint of W. 
We recall that the conditional densities of SRB measures are Holder continuous, and their 
fluctuations are uniformly bounded [5l [TO] . Now eq. ([38]) implies 

m{R e>s ) < C"5 

for some constant C" > 0, hence p{9) is uniformly bounded. 
Next we prove that p{9) is continuous. For 9\ ~ 9 we have 

p{9 x ) - p(9) = Inn 

Now there are unstable manifolds that cross both rectangles Re lt s and Rgj and those which 
cross only one of them; accordingly we have 

p(Re u5 ) - p(Re,s) = Ai + A 2 , 

where Ai accounts for the former, and A2 for the latter. The conditional density of p on 
each unstable manifold is Holder continuous, and unstable manifolds have uniformly bounded 
curvature [U [TO] , hence once can easily see that 

I Ai| < C<5|0i -#| 7 (39) 
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for some constants C > and 7 > (in fact, 7 = 1/3 for our model; cf. [5j Corollary 5.30]). 
It remains to show that 

lim lim sup | A 2 1 18 = , (40) 

i.e., the contribution from unstable manifolds crossing just one rectangle is negligible. 

To estimate A2, denote by F(y) = p(UW: \W\ < y) be the measure of all the unstable 
manifolds of length < y. Then eq. (|38H can be written as 



F{28) + / — dF{y) < CS. 

J25 y 



Dividing by 5 and taking the limit 5 — > gives 



00 



2 1 ^1<C (41) 



'0 y 

(see also [5j Exercise 7.15]). Now let Lg = {(6,ip): tp G [— tt/2, vr/2]} denote the line in M. 
with the fixed 9 coordinate. Denote 

Fe,e 1 (y) = p(L)W: \W\ < y, W terminates between Lg and Lg^. 

Then we have 

to^<fM (42) 

5^0 ° Jo y 

Now, as 9\ and 6 get closer together, Fg^g^y) monotonically decreases for each fixed y > 0. 
Moreover, we have 

\imF m (y) = Vy>0 (43) 

01-S-0 

because the union of unstable manifolds terminating exactly on the line Lg has //-measure 
zero. To see this observe that unstable manifolds terminate on singularity lines of the past 
iterations of the collision map T, i.e., on singularity lines of T~ n , n > 0. These lines inter- 
sect the line Lg at countably many points thus there are at most countably many unstable 
manifolds terminating on Lg. Finally each individual unstable manifold has fi- measure zero, 
as is guaranteed by the Poincare recurrence theorem. Now combining eq. (l4ip -eq.( l43p proves 
eq.(g0D. 

In smooth hyperbolic systems without singularities all unstable manifolds are long enough 
so that A2 = 0. Then the density p is not only continuous, but Holder continuous, according 
to eq. (|39p . We believe that in our MH model, too, the main contribution to the structure of 
p comes from Ai, so that p is also Holder continuous, but our estimate eq. (l40p on A2 is too 
poor to prove that. 

The above argument applies to projections of SRB measures onto some coordinates 
(transversal to stable and unstable directions). If, instead of projections, we integrate smooth 
functions like in eq. (|14p . then those functions can be incorporated into conditional densities 
of the SRB measure on unstable manifolds, and the argument will work for that situation, 
too. 



B Derivation of Green-Kubo formulas 

The derivation of eqs.([9l fl~0|) and (|18p is based on a Kawasaki- type formula used in linear 
response theory, see [4]. For a small external field E and the corresponding SRB measure 
/tg, we can integrate any bounded piecewise Holder continuous function Je on M. as follows: 

00 

HeUe) = Me) ~ e ~ EAE ))> (44) 

n=l 



18 



where Ae denotes the displacement of the particle in the direction of the field (i.e., in the 
positive x direction) during its free flight to the next collision. More precisely, for (9, ijj) € M 
we set 

A E {9^) = / cos($f(M))dt (45) 

J o 

where we use this definition to avoid the ambiguity on the difference between two points on 
a torus (periodic boundary conditions). The above Kawasaki- type formula eq. (|44p is derived 
in [I Eq. (16)]. 

The function Je may depend on the field E, but it must have a limit /o = lim/g as 
E — > 0. Since To = lim^o Te and to = lime^o te, as well as Ao = lime_j.o Ag = r cos (f>, a 
first order Taylor expansion of the infinite sum in eq. (!44p gives 

oo 

HeVe) = MJe) +i?^ W ((/o o T n )A ) + o(E), (46) 

n=l 

see [H Eq. (17)]. The infinite sum in the above equation converges because //o((/o°T n )Ao) — > 
A*o(/o)/ Lt o(Ao) = exponentially fast as n — > oo for fo Holder continuous (decay of correla- 
tions) and //o(Ao) = 0. Note that this sum is independent of E, hence the second term is 
linear in E. The first term /io(/fi) will be handled separately for each Je- 

B.l Derivation of Eqs.flS) and flTD]) 

We first derive eq.© for the local density rag(x), which can be represented, according to 
eq-©, by 

<5 2 n s (x) = MJa) + o(S 2 ) (47) 

where A = {x': ||x' — x||oo < 6/2} denotes the square in Q with side 5 centered on x. We fix 
a small 5 > and will estimate the integral in eq. (|47|) . 

It is convenient to extend the space M by adding the sides of the square A to dQ. In other 
words, every time a trajectory crosses the boundary of A and enters or exits A we register a 
'virtual collision' (the trajectory does not actually change its direction, so the collision has no 
effect on the trajectory of the particle, i.e. dA plays the role of 'transparent walls'). Adding 
transparent walls with virtual collisions is a useful trick in the study of billiards. 

By adding this transparent wall the phase space of the system is thus extended from A4 = 
[0,47r] x [-7r/2,7r/2] to Ma = [0,47r + 4<5] x [-tt/2,tt/2] where 6 G [4tt, 4tt + 4<5] parametrize 
dA. Consistently we must replace the map Te with the new map Te s a constructed as in 
section [1] the SRB measure /ig with fiE,A defined as in eq.© and the function with 
J®' defined as in eq.©. The notation jj^' helps keeping track of the fact the A appears 
both in the indicator function appearing in eq.© and in the phase space Ma on which J^' A 
is defined. From this follows that J A ' ^ if and only if 9 represents a collision taking 
place on dA and the outgoing velocity points inside A. Indeed, /^(X^fl, tp)) ^ if and only 
if X((^,^) G A. But in this case the last collision of the trajectory was with dA and the 
velocity was pointing inside A. 

Clearly eq. ([4*6j) remains true if we replace he with he, A- We apply it to /e = J a' an d 

get 

oo 

5 2 n E {*) = Ho,a(Ja A ) + I^Axa) + E £ ho,a{( j a A t o,a) A o,a) + o{E) (48) 

n=l 

where Ao,a is defined as in eq. (f45j) and Xa = ^a' A ~ ^A. A - direct calculations we get 

^ (J ^ ) = Ar^(Q) =52no(x) - 
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Observe that, for every E, J A ' A = O(S) while //o,A(supp( J A ,A )) = 0(8). Moreover, from 
eq.© we have that |Xf (6,ip) - X?(0,^)| = 0(Et 2 ) so that J A ' A - 4' A = 0(5 2 £). Thus 
even though the term jUo,a(x]|) = 0(E8 3 ) is linear in i?, its contribution vanishes in the limit 
5—t-O and we will ignore it. We finally arrive at 

oo 

n s (x) = n (x) + 8~ 2 E ]T ^A^Y V>,a)*o) + °(^)- ( 49 ) 

n=l 

Observe that a trajectory originating from a point x gives a non zero contribution to the n-th 
term in the sum only if its n-th collision is with dA and the trajectory enters A. Since E and 
8 are small, this implies that the n — 1-th collision was with a non-virtual obstacle. We can 
thus rewrite eq. (|49p as 

d(x) = <T 2 L Q (J A A 0>A ) +f>oM°T ")A )J (50) 

where we have neglected the trajectories that collide more than once with A, since they 
contribute 0(8 3 ) to the integral, and the difference between Te and te,a, appearing in J A ' A , 
since it is 0(E) and thus does not contribute at first order. 

We still have to discuss the limit 8 — > in eq. (|50p . This limit is non trivial since, although 
the correlations appearing in the infinite sum decay exponentially for every A, we need to 
show that such a decay is uniform in 8. To show this we can take the limit <5 — > term by 
term in the sum. Taking this into account we obtain eq.©. 

We can now derive eq. (fT0|) for the local average velocity vg(x). Similarly to eq. (|4"T|) we 
have 

8 2 ^) = ^^E(Hf) + o(8 2 ), (51) 

where 

Hf (0,</>) = — / (cos($f (9^)),sm(^f (6^)))I A (X t E (9,^)) dt (52) 

T E Jo 

and A again denotes the square in Q with side 8 centered on x. In this case we will not need 
to introduce virtual collisions with A like we did for eq.(|9]). Applying eq. (|44p we obtain 

oo 

<5 2 n E (x)v fi (x) = W (Hf) + £^/i ((H° oT")A ) + o(E) (53) 

n=l 

To eliminate the first term //q(H^) we apply an antisymmetrization. Due to the invariance 
of he we have /j,e(H a ) = ue(H^ o T^ 1 ), hence, by applying eq. (j4"4"j) to H A o Tg , we get 

oo 

<5 2 n s (x)v £ (x) = f i (UioT- 1 ) + EY,^{(H.A°n i )A ) + o(E) (54) 

n=0 

Next, let J: M —> M denote an involution defined by J(0,tfj) = C(9,—tp). Due to the 
time reversibility of the perturbed dynamics we have J o Te = T^ 1 o J , and therefore 
H A = — H A o J o Te- Also note that fig is invariant under both To and J . Thus if we add 
eq. lfollj) and eq.(|54p. the first terms cancel out. 

Moreover, the time reversibility of the billiard dynamics implies Ao = — Ao o J o To and 
H A o T n = -H^ o T ~ n o J o T for all n. Therefore, 

(H° A o T n ) • A = [(H^ o V*) A ] o J o T . 
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Thus adding eq. (|57l) and eq. (|58j) together gives 

oo 
■n.= — oo 

Taking into account eq.(@9]) we get 

oo 

M*) = 2F^W E Vo((K°a o T «)A ) (55) 

7I = — OG 

where, due to the time reversibility of the dynamics, the term for n and — n in the sum are 
equal. To obtain eq. (fT0j) we have used that lim^o ^~ 2 Mo((H^ oT^)A ) = c(p£(A ), p s n (A )) 
with c = no(x). 



B.2 Derivation of Eq.flTS) 



To derive eq. (fT8|) , we apply eq. PB"j) to /e = J|v^ ^ , which was defined in eq. ([6]) and where 
A(4>,5) is the set <f> — 5/2 < 6 + ip < <p + 5/2, i.e. the set of velocity vectors that make an 
angle <f>' G — 5/2, eft + 5/2]. Note that when = 0, the trajectory is a straight line, so 
/o = 70-^4(0,5) • 

As in eq.flT]) we can write 5h E {(t)) = /"e(-^a(0,5)) + °(^)- We denote A^ = A • lA(<j>,6) 
and recall that Ao = tocos</>, so that Jams) COS( l ) = A^/fo. Thus eq. (fl6j) becomes 

E °° 

5 cos <j> h E {4>) = w{Ja{4>,5)) cos <t> + — E M( A <M o T n )A ) (56) 

^ n=l 

and, again due to the invariance of he, we get 



oo 



E 

5 cos <t>h E {<t>) = Mo(^f(0,(i) °7^ 1 )cos(/>+ — ^^ ((A^(5 o T n )A ) (57) 

T ° n=0 

and (denoting, for brevity, 0~ = (f> + n) 



oo 



Scoscf)- h E {<jT) = /i (Jf^- 5) )cos<T + -J2^{(^-,s°U l )A ), (58) 

" n=l 

In all the above formulas we have suppressed the o{E5) terms. Time reversibility implies 
^A((f>5) -^B 1 = "^4(0- <5) ° wnere <^ denotes the time reversal involution. Thus, if we add 
eq. (|57p and eq. (|58|) . their first terms cancel out. Similarly, we get Ao = — Ao ° J o Tq and 
A<a-,(5 ° lo 1 = — A0 j( 5 o T ~ n o J o Tq for all n. Therefore, 

(A - , o T ") • A = [(A^ o T~ n ) A ] oJoT . 

Thus adding eq. ([57|) and eq. (|55j) gives 

M£>-M0 + 7r) = - t E M( A ^s°T n )A )+o(E). (59) 

OTq COS ^-^ 



n=— oo 



This is an infinite sum of correlations which decay exponentially fast [5]. 
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In the case of the billiard shown in Figure [T] we have that, due to the symmetry of the 
system, + 7r) = so that we get 

^ oo 

a ^ = vF 2 E Mo((A^oT ")A )) 

ZOTq cos m 

n=—oa 
1 oo 

= ^F 2 E w(A^-(A oT ")) 

iOTn COS <p 

n=—oo 

where we used the invariance of (jlq under To. Finally using the relation A^ = tqI^ ^ cos0 
and taking the limit 5 — > gives eq. (|19p . 



C Isotropy of the Collision Time 

As a motivation for our definition of /j,q(tq), we define a directional mean free time as follows. 
Given 4> £ [0, 2tt] and 5 > 0, let / = Iau,S) be as in Section IB. 21 Due to ergodicity, we have 

lim Er=oVoCT,V0)/(r^)) = MtqJ) 

for almost every (0,ip) € .M, and we call the limit (if it exists) 

= lim — — 

the directional mean free time (corresponding to the angle 0). Now arguing as in Section TB.2I 
we get 

lim <TVo(toI) = I t{6- 6, 6) cos(0 - 0)x(9) d8 (61) 
<5^o J 

and 

lim 5~Vo(^) = / cos((P - 6)x{0) d6 = Z+, (62) 

recall eq. (|17p . Therefore 



^ = /Wo) 

is the conditional expectation of r defined in Section 12.21 

Now it is easy to see that cos(</>— 9)x(d) d6 is the length element in the direction orthogonal 
to the outgoing velocity vector (i.e., in the direction </>+7r/2). Therefore the integral in eq. (|61|) 
is equal to the area of the billiard table, which is 1 — it{r\ + r|) in our case. Thus 

6 1 - 7r(rf + r\) 

2(n + r 2 ) 

which is constant (independent of 4>). 

The above argument generalizes to any Sinai billiard with convex obstacles B±, . . . ,B p . 
For each obstacle Bk and angle <f> we denote by width^(£?fc) the "width" of Bk in the direction 
orthogonal to <f), i.e., the length of the projection of B^ onto a line orthogonal to all velocities 
running at the angle <j). Then by the above argument we have 

Area(Q) 

T * E.width^r [M) 
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This formula holds for each </> except directions in which billiard trajectories can run indefi- 
nitely without collisions. If the horizon is finite, no such trajectory exists, and eq. ()63p holds 
for every <p. If the obstacles are circular disks, as they are in our studies, the "width" of 
.Bfc is just its diameter, and comparing eq. ([63|) with eq. ([64"j) we see that is constant, i.e., 
independent of </>. 

It is not hard to see that averaging fU over (ft gives 



1 

2^ 



T^dtf) = t = hq{t ), 



the classical (unconditional) mean free path, which is known to be 



- _ 7T ■ Area(Q) 

length(dQ) ' 1 ' 



see (5j Section 2.13]. If is constant, then of course = tq for all 
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